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Abstract. We study the classical and quantum perturbation theory for two non-resonant oscillators 
coupled by a nonlinear quartic interaction. In particular we analyze the question of quantum corrections 
to the torus quantization of the classical perturbation theory (semiclassical mechanics) . We obtain up 
to the second order of perturbation theory an explicit analytical formula for the quantum energy levels, 
which is the semiclassical one plus quantum corrections. We compare the "exact" quantum levels obtained 
numerically to the semiclassical levels studying also the effects of quantum corrections. 
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1 Introduction 

Nowadays there is considerable renewed interest in the transition from classical mechanics to quantum 
mechanics, a powerful motivation behind that being the problem of the so-called quantum chaos [1-3]. 
An important aspect is represented by the semiclassical quantization formula of the (regular) energy 
levels for quasi-integrable systems [4-6], the so-called torus quantization, initiated by Einstein [7] and 
completed by Maslov [8] . 

It has been recently shown [9,10] that, for perturbed non-resonant harmonic oscillators, the algorithm 
of classical perturbation theory can be used to formulate the quantum mechanical perturbation theory 
as the semiclassically quantized classical perturbation theory equipped with the quantum corrections in 
powers of fi "correcting" the classical Hamiltonian that appears in the classical algorithm. In effect, one 
can explicitly calculate corrections to the Einstein-Brillouin-Keller (EBK) quantization of the classical 
tori [4]. For example, the quantum corrections of the one-dimensional x 4 perturbed harmonic oscillator 
have been studied in great details by Alvarez, Gram and Silverstone [11]. 

Examples of rather detailed studies of semiclassical approximations, their resummations and of the 
Birkkoff-Gustavson normal forms can be found in Ali and Wood [12] and Ali, Wood and Devitt [13]. 
Another example of a rather complete semiclassical analysis of a one-dimensional system, namely the 
quartic oscillator, has been published by Voros [14]. 

The aim of this paper is to extend previous studies to a two-dimensional system, which is more 
interesting because it is nonintegrable and thus generic. The integrable systems are rather exceptional in 
the sense that they are typically isolated points in the functional space of Hamiltonians and their measure 
is zero in this space. If we randomly choose a system realized in nature, the probability is one that the 
system is nonintegrable [15]. 



2 Classical Perturbation Theory 

The model is given by two non-resonant oscillators coupled by a nonlinear quartic interaction of strength 
g (Pullen and Edmonds [16]): 

^ = y(P? + ??) + y(^ + ^)+54 2 (1) 



Note that the same Hamiltonian has been obtained for the Yang-Mills-Higgs classical mechanics. The 
non-intcgrability of this system and its transition from order to chaos has been studied in classical with 
Poincare sections and in quantum mechanics with the spacing distribution of energy levels [17,18]. 
Through the canonical transformation in action-angle variables [19,20]: 

q k = V24 cos6» fc , Pk = V^ifesinflfe, k = 1,2 (2) 
the Hamiltonian can be written: 

H = H o (I 1 ,I 2 )+gV(h,I 2 ,0 1 ,0 2 ), (3) 

where: 

H (h,I 2 ) =w 1 h+w 2 I 2 , (4) 

V(I U 12,0!, 2 ) = 4/i/ 2 COS 2 9 l COS 2 e 2 . (5) 

Following the classical perturbation theory [19,20], we search for a canonical transformation (h,I 2 , 0\,0 2 ) 
(h,I 2 ,6i,6 2 ) to obtain a new Hamiltonian that depends only on the new action variables up to the second 
order in a power series of g: 

H{h,I 2 ) = H {h,I 2 ) + gHi(hJ 2 ) + g 2 H 2 (h,I 2 ). (6) 

The generator S of the canonical transformation is supposed to be capable of expansion as a power 
series in g of the form: 

s(h ,i 2 ,8i, 02) = hex + i 2 o 2 + 9 s 1 (h ,12,0!, e 2 ) + g 2 s 2 (h , i 2 , e u e 2 ), (7) 

and to satisfy the equations: 

dS ~ dSi 2 8S 2 

Ik = w k =Ik+g w k +9 w k > (8) 

2 dS a dSi 2 dS 2 

°k = 1 r r = k + g- r +g 2 - r . (9) 
ol k dl k dl k 

From the Hamilton-Jacobi equation: 

H o(^ W2 )+9V( Wi , W2 ,e 1 ,9 2 ) = H a (I 1 ,I 2 )+gH 1 (I 1 ,I 2 )+g 2 H 2 (I 1 ,I 2 ), (10) 

we have a number of differential equations that result on equating the coefficients of the powers of g: 

H (h,I 2 ) = H (h,I 2 )=Lu 1 h+uj 2 I 2 , (11) 

H 1 (I 1 J 2 ) = (u> 1 ^+W2-^) + V{I 1 ,h,e 1 ,02), (12) 

u(f t\ ( ds ^ ds 2 dv ds 1 dv os, 

H 2 (I U I 2 ) = {U1 — +Ua —) + (— — + — — ). (13) 

The unknown functions Hi, S\, H 2 and 52 may be determined by averaging over the time variation of 

the unperturbed motion. 

At the first order in g we obtain: 

J 2 ) - 7^2 J o J o d9 1 d0 2 V{h,I 2 ,0 U 2 ) = hh, (14) 

and 

S 1 (I 1 ,I 2 ,0 1 ,e 2 ) = -7/1/2 [—sin 20! + —sm20 2 + 

4 UOi U! 2 

+ sin 2(0! - 2 ) + ; sin 2{9 X + 9 2 )}. (15) 

At the second order in g we have: 

H ^ h) = ^L Jo ^^dhW.+dhW^ 



. jM li{ k + k ) .(J^M + ik±M ] . (16) 

8 UJ 2 U>1 Ml — UJ 2 OJl+ UJ 2 

We observe that the integrable approximate Hamiltonian (6) could be obtained alternatively as the 
Birkhoff-Gustavson normal form, which is a purely algebraic method of calculating the action variables 
order by order for perturbed harmonic oscillators with polynomial perturbations [21]. The Hamiltonian 
(6) depends only on the actions so that by an application of the EBK rule: 

h = (ni + ^)h, h = (n 2 + -)h, (17) 

to equations (11), (14) and (16), we obtain a semiclassical analytical formula of the energy levels. Obvi- 
ously h is the Planck constant and n\, n 2 are two integer quantum number. 

Note that the torus quantization can be applied because the normal form obtained (equation 6) is 
integrable. The new action variables remain very close to their initial conditions. This is no longer true, in 
general, for higher dimensions (see [19,20]). In the next section we show how to connect our semiclassical 
formula with the usual quantum perturbation theory. 



3 Quantum Perturbation Theory 

In quantum mechanics the generalized coordinates are operators which satisfy the usual commutation 
rules for Bosons [4k, Pi] = q k pi — Pi4k = indki, with k,l = 1,2. Introducing the creation and destruction 
operators: 

afe = -7= (4k +iPk), °i = —={q k - ip k ), (18) 
V in V 2n 



the quantum Hamiltonian can be written: 

H = H +gV, (19) 

where: 

H = fuv 1 (afai + -) + Jkv 2 (a£a 2 + -), (20) 

V = ^-(a 1+ a+) 2 (a 2 +a+) 2 . (21) 
If \n\n 2 > is the basis of the occupation numbers of the two harmonic oscillators, the matrix elements 

are: 

< n[n 2 \H \n 1 n 2 >= %i(m + -)+u 2 (m + ^R'^A^ ( 22 ) 

and: 

< ni n 2 \V\ ni n 2 >= — [V»i(ni - l)J„ ini _ 2 + y/(m + l)(m + 2)S n[ni+2 + (2m + l)<5„ ; „Jx 

x Wn 2 (n 2 - l)S <n2 _ 2 + V(n 2 + l)(n 2 + 2)S <n2+2 + (2n 2 + l)5 <na ]. (23) 
The Rayleigh-Schrodingcr perturbation theory [22] up to the second order gives us: 

E(mh, n 2 h) = E (mh, n 2 h) + gE^mfr, n 2 h) + g 2 E 2 (n 1 h, n 2 h), (24) 



where: 



E Q {n 1 Ti,n 2 h) = h[ui(m + i) +u 2 {n 1 + i)], (25) 



Ei(mh,n 2 h) —< nin 2 \V\nin 2 >, (26) 

| < n^^jVjning > | 2 _ 
K[u>i(ni - n'j) + u> 2 (n 2 - n' 2 )] 



E ( h h) ^ < n i n 2l^l n i n 2 > | 2 ^27~j 



(n' 1 ,rt2)7^(ni,n 2 ) 



We obtain immediately: 



and after some calculations: 



E 1 {n 1 h,n 2 h)=h 2 {n 1 + ^)(n 2 + ^), (28) 



„ , n n ^ fi 3 r ni(TH-l)n 2 (Ti 2 -l) (n 1 + l)(n 1 + 2)(n 2 + l)(n 2 + 2) 

ti 2 (nih,n 2 h) = — ■ 1 h 

32 cji + w 2 Wi + w 2 

| m(ni - l)(ra 2 + l)(rt 2 + 2) _ (rti + l)(ni + 2)n 2 (n 2 - 1) 
Wi — w 2 Wi — w 2 

| ni(m - l)(2n 2 + l) 2 _ (ni + l)(m +2)(2n 2 + l) 2 

(2m + l) 2 n 2 (n 2 - 1) _ (2m + l) 2 (n 2 + l)(n 2 + 2) g 

The zero and first order quantum terms coincide with the semi-classical ones: 

E ( ni n, n 2 Ti) = H {{ni + ^)h, (n 2 + (30) 

Etim^mh) = Hi((m + ±)h, {n 2 + (31) 

and the second order quantum term can be written: 

E 2 ( ni n, n 2 K) = H 2 {{m + ^)h, (n 2 + + h 2 Q 2 {{m + (n 2 + (32) 

where: 

n » . 1^ , . I-.*-. 3 r (m-n 2 )ft (ni+n 2 + l)^ 

Q 2 ((m + x)?i, (n 2 + = - — + — — ; ]■ (33) 

Z Z 6Z LOl — U) 2 U)\ "T U) 2 

The quantum series rearranges directly into the classical canonical perturbation series plus quantum 
corrections proportional to successive powers of h [10]. 

The term h 2 Q 2 represents the quantum corrections to the EBK quantization up to the second order 
of perturbation theory. These quantum corrections depend linearly on quantum numbers. For a quite 
similar system, Robnik [21] obtained that the spectra differ only by an additive constant independent of 
the quantum numbers (actions), but he studied only the first order of perturbation theory. 

To conclude this section we mention the general problem of quantization: Only the quantization of co- 
ordinate space or the quantization of linear canonical transformations of coordinate space yields the right 
quantum mechanics, whose results agree with experiments. This is exactly the approach implemented in 
our present case (18), and thus our quantization is equivalent to the coordinate space quantization. 



4 Numerical Results 

We compute the energy levels with a numerical diagonalization of the truncated matrix of the Hamiltonian 
(19) in the basis of the unperturbed oscillators (see [23] for a more technical discussion). The numerical 
energy levels depend on the dimension of the truncated matrix: We compute the numerical levels in 
double precision increasing the matrix dimension until the first 100 levels converge within 8 digits (matrix 
dimension 1225 x 1225). This means that the agreement between numerical values and real ones is very 
good and we use the word "exact" to indicate the numerically computed energy levels. 

In Table 1 we compare the "exact" levels to the levels of the quantum perturbation theory and to the 
semiclassical results. A very good agreement is observed for the lowest energy levels. 

In table 2 we show the error in units of the mean level spacing D between the "exact" levels and the 
levels obtained with semiclassical and quantum perturbation theory. We observe that the algorithm pro- 
vided by the appropriate semiclassical quantization is comparable to the algorithm provided by ordinary 
quantum perturbation theory but the quantal corrections do not always increase the accuracy. Thus 
for some of the calculated levels the semiclassical quantization gives better results than the quantum 
perturbation theory. 



In Table 3 we show the error, in units of the mean level spacing D, between the "exact" levels and 
the semiclassical levels. As is seen, by decreasing T% from 1 to 1CP 1 the quality of the approximation 
improves considerably. However, it is important to observe that if Ti, no matter how small, is kept fixed, 
the EBK quantization (torus quantization) of the individual levels is only a first order approximation of 
an expansion in Ti. Therefore, in general the accuracy of the approximation decreases for higher levels. 
To get a good agreement it is necessary, as is well known, to implement the classical limit, i.e. % — > and 
m, ri2 — > oo, while at the same time keeping the actions I\ = (m + 1/2)71. and I 2 = (n 2 + l/2)h constant 
[23]. 

5 Conclusions 

We have examined the transition between the classical and the quantum mechanics for a two-dimensional, 
nonintegrable and non-resonant system. Up to the second order of perturbation theory we have decom- 
posed the quantum description into the classical description (i.e. the leading semiclassical term) plus 
quantum corrections which depend linearly on the quantum numbers. The semiclassical energy levels 
and the levels obtained with quantum perturbation theory are in good agreement with the "exact" nu- 
merical ones: The semiclassical quantization is comparable to the quantum perturbation theory, and for 
some levels the semiclassical quantization gives better results than quantum perturbation theory. 

There arc two series involved in the work. One in powers of K because of the semiclassical approxi- 
mation and another in the interaction parameter g. The first one is an asymptotic expansion (see some 
very recent results in [24-27]) but the second one may be not convergent [28]. The classical and quantum 
perturbation series typically diverge and thus do not necessarily describe the exact levels, not even after 
a certain resummation (except for some important notable exceptions like the anharmonic oscillators 
with / degrees of freedom with a polynomial perturbing potential which is asymptotically positive defi- 
nite [29,30]), it is important to compare the semiclassical approximation (and the quantal perturbation 
results) with the exact spectra, which in general is impossible, since we generally do not have explicit 
solutions of the Schrodingcr problem in a closed form. Therefore, we stress the importance of specific 
case studies like the present one, in order to get a better understanding of the quality of semiclassical 
mechanics. 

Finally we note that the extraction of quantum corrections for resonant systems is a more intricate 
procedure; some initial results for perturbed resonant oscillators can be found in Graffi [31]. 
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055694 


7 
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389530 
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093505 
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088912 


(2,3) 



Table 1: Comparison between "exact" levels and levels obtained by perturbation theories. First 20 
levels. E ex are "exact" levels, E^ c iJl2 are semiclassical levels, and £^ nj are levels obtained with quantum 
perturbation theory, where n\ and n 2 are the quantum number. % = 1, g = 10 _1 , uj\ = 1 and uj 2 = V%- 
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(2,0) 


3.6215890-10- 2 


4.2791300-10 
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Table 2: The error measured in units of the mean level spacing D for the first 20 levels. D is 
calculated for the lowest 100 levels. E ex are "exact" levels, £^„ 2 are semiclassical levels, and £^„ 2 are 
levels obtained with quantum perturbation theory, where n\ and n 2 are the quantum number. % = 1, 
g = 10 _1 , u)\ = 1 and uu 2 = V%- 
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Table 3: The error measured in units of the mean level spacing D between "exact" levels and 
semiclassical levels. First 20 levels. D is calculated for the lowest 100 levels, g = 10 _1 , u\ = 1 and 

LU2 = V2. 
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